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Abstract. A recent experiment showed that cylindrical segments of water filling a 
hydrophilic stripe on an otherwise hydrophobic surface display a capillary instability 
when their volume is increased beyond the critical volume at which their apparent 
contact angle on the surface reaches ninety degrees (Gau et ai, Science, 283, 1999). 
Surprisingly, the fluid segments did not break up into droplets — as would be expected 
for a classical Rayleigh-Plateau instability — but instead displayed a long-wavelength 
instability where all excess fluid gathered in a single bulge along each stripe. We 
consider here the dynamics of the flow instability associated with this setup. We 
perform a linear stability analysis of the capillary flow problem in the inviscid limit. 
We first confirm previous work showing that that all cylindrical segments are linearly 
unstable if (and only if) their apparent contact angle is larger than ninety degrees. 
We then demonstrate that the most unstable wavenumber for the surface perturbation 
decreases to zero as the apparent contact angle of the fluid on the surface approaches 
ninety degrees, allowing us to re-interpret the creation of bulges in the experiment as 
a zero-wavenumber capillary instability. A variation of the stability calculation is also 
considered for the case of a hydrophilic stripe located on a wedge-like geometry. 
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Capillary instabilities are phenomena we witness in our daily lives, and their study is a 
field with a rich history [l|, y, 0, 0] ■ The classical Rayleigh-Plateau instability refers to the 
surface-tension induced instability of a cylindrical liquid column. For volume-preserving 
deformations of sufficiently long wavelengths along a fiuid cylinder, the surface area of 
the fiuid can be made to decrease. These deformations lower the surface energy of the 
fluid and are therefore favorable, so an infinite cylindrical fiuid column is always capillary 
unstable. For a cylindrical fiuid column of radius R, density p and surface tension 7, 
this instability is of inviscid nature, and occurs on a typical time scale ti ~ {pRV'yY^'^, 
with the most unstable wavelength being on the order of the column radius 0] • 

Many variations on this classical result have been considered in the past, and 
we refer to Ref. jl!] for a review. In the present paper, we consider such instabilities 
when they occur for a cylindrical segment of fiuid in contact with a solid. A recent 
experimental investigation on the stability of a cylindrical segment of liquid on a stripe 
of hydrophilic material on an otherwise hydrophobic surface has shown surprising long- 
wavelength instabilities. In that case, the fiuid segment was observed not to break up 
in many droplets as in the classical Rayleigh-Plateau instability, but instead to form 
a single large bulge [51 . The purpose of this paper is to study the dynamics of such a 
surprising capillary instability. 

A number of previous studies have approached the problem of fiow stability for a 
cylindrical segment of fiuid pinned on a solid substrate. Davis considered a fiow rivulet 
down an incline, with its contact line pinned along a stripe jsl . Using an energy method, 
it was found that when the interior angle {i.e. the apparent contact angle of the rivulet 
on the horizontal surface) was less than 90°, the rivulet is stable, and it is unstable 
otherwise. A similar result was later obtained using energy minimization considerations 




(a) (b) 

Figure 1. The growing instability of a cylindrical segment on a hydrophilic stripe. 
Left: For small enough volume, the cylindrical segments, located along the hydrophilic 
stripe and separated by hydrophobic stripes, are stable. Right: At a critical value of 
the fluid volume, corresponding to an apparent contact angle of the fluid of ninety 
degrees, the cylindrical segment become unstable, and solitary bulges grow. Reprinted 
from Gau et al, Science, 283, 1999 3, courtesy of R. Lipowsky and with permission 
from AAAS doi : 10 . 1126/ science . 283 . 53 98 . 46 . 
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(a) (b) (c) (d) (e) 

Figure 2. Schematic representation of the different stages of growth of a cylindrical 
segment on a hydrophilic stripe in the absence of any capillary instability, (a) For 
small enough volumes, the fluid grows solely on the hydrophilic region, and displays 
an advancing contact angle, 0^; (b) At a critical fluid volume, the line separating 
hydrophilic and hydrophobic domains in reached; (c) As the fluid volume continues to 
increase, the contact line remains pinned, and the apparent angle increases; (d) The 
angle stops increasing when the apparent angle reaches the value, 9'^, of the advancing 
angle on the hydrophobic domain (note that 9^ > 90° > 9\)] (e) As the fluid volume 
continues to increase, the contact line region now moves into the hydrophobic domains, 
with angle 9^. In this paper, we show that the cylindrical segment displays capillary 
instability as soon as (c) is reached (apparent angle of 90°). 



for a cylindrical interface pinned to a slot [71]. In addition, the critical length for the 
instability of a finite cylinder was found to become infinite at the critical angle of 
ninety degrees [7|. Similar results were later recovered using differential geometry, 
where volume-preserving perturbations were seen to lead to decrease in surface area 
of the cylindrical filament only in the case where the apparent angle exceeded 90°, and 
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The experimental observations of Gau et al. [5[, which motivate the present study, 
are reproduced in Fig. [H The fluid (water) is filling hydrophilic stripes (contact angle 
~ 5°) on an otherwise hydrophobic substrate (contact angle ~ 108°). The fluid volume 
is then increased. In the absence of the observed instability, a schematic representation 
of the growth process is represented in Fig. [2j As the fluid volume increases, the contact 
line at the edges of the cylindrical segments of water remains pinned, until the volume 
of fluid becomes large enough that the apparent contact angle is equal to the advancing 
contact angle on the hydrophobic surface. At this point, any subsequent volume change 
would be accompanied by a motion of the contact line into the hydrophobic substrate. 
The experimental result obtained by Gau et al. is that such a process is unstable. As 
soon at the apparent contact angle of the fluid segment reaches the critical value of 
90° (intermediate between the contact angles on both surfaces; see Fig. [2]3), the surface 
becomes capillary unstable. This result is consistent with the previous studies discussed 



becomes capillarv 
above 0, H H, 0, 



10, lUj. Surprisingly, and in contrast with a free fluid cylinder, the 
observed unstable mode does not display a wavelength on the order of the cyhndrical 
radius, but instead the wavelength appears to be much larger, and the segment evolves 
to a state where all the excess fluid gathers on single bulge (Fig. [1]). Some features of 
these experiments were later reproduced by Darhuber et al. in their experimental and 
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consistent with the simplified one- dimensional stability study of Ref. 

The aim of the present work is to focus on the dynamics of the instability process 
by performing a linear stability analysis of the cylindrical segment in the experiments 
of Gau et aL, and predicting the dependence of the growth rates and most unstable 
wavelengths of the unstable modes on the apparent contact angle of the cylindrical 
segment. We restrict our study to low viscosity liquids, such as water, and perform an 
inviscid study. Physically, inviscid capillary instabilities of a fluid column are expected 
to occur on a time scale ri ~ {pR^ /'jY^'^. For comparison, the time scale for viscous 
effects to propagate diffusively across the width of a fluid column is T2 ~ R"^/ v where v is 
the kinematic viscosity of the fluid. The inviscid approach will therefore be a reasonable 
modeling assumption as long as Ti <^ r2, which is equivalent to i? ^ with i = pz/^/o" 
is the Ohnesorge length scale of the fluid. For water, i is on the tens of nanometers, 
which is much smaller than the typical cross sectional size in the experiments of Gau et 
al. jsf (tens of microns; see Fig. [1]). As we show below, within these assumptions, the 
most unstable wavelength for the capillary instability of the cylindrical segment tends 
toward infinity as the contact angle approaches 90°, thereby allowing us to re-interpret 
the creation of bulges in the experiment as a zero-wavenumber capillary instability [sl . 

2. Setup and linear stability 

The geometrical setup for our linear stability calculation is illustrated in Fig. [31 The 
basic state is a cylindrical segment of fluid of radius R, whose two-dimensional contact 
line is pinned along a stripe. The apparent contact angle of the fluid at the contact line 
is denoted 6c. We assume dc to sufficiently larger (smaller) than the advancing angle on 
the hydrophihc (hydrophobic) substrate so that we can safely assume that during the 
initial stages of the instability the contact line remains pinned along the same location. 
In addition, since we know from earlier work that the case where 6c < 90° is stable 
we will focus here on determining the dynamics of the instability in 



the case where 9r > 90°. 



2.1. Governing equations and linearization 

Assuming the liquid is both incompressible and inviscid, its motion is prescribed by the 
Euler equation and the continuity equation, 

p(^ + u-Vu] = -Vp, (1) 



.dt 

V ■ u = 0, (2) 

where u and p denote the velocity and pressure fields respectively. Since there is no flow 
in the basic state, the linearized equations are simply 



Capillary instability on a hydrophilic stripe 



5 




Figure 3. Setup for the inviscid stability calculation. We consider a cylindrical 
segment with a pinned contact line, displaying an apparent contact angle 6c > 90°. 
The radius of the cylindrical segment is denoted R, and the height of its center about 
the surface is denoted h. 



V ■ u' = 0, (4) 

where primes denotes small deviations from the basic state. We employ cylindrical 
coordinates, and use the center of the cylindrical segment as the origin of our coordinate 
system; the variable z denotes therefore the coordinate along the stripe, and 9 = 
denotes in the vertical direction. Let the radius of the free surface be parameterized as 
C,{t, 6, z). Since the radius of the free surface is R in the basic state, then the perturbation 
in the position of the free surface is C,' = C, — R- The free surface moves with the local 
velocity of the fluid, which after linearization is written as 

<\r^R = ^ • (5) 

Additionally, at the contact points, the position of the contact point is fixed, giving the 
condition 

Ur\{r=R,e=±ec} = 0- (6) 
The jump in pressure across the free surface is related to the curvature of the surface 
and the surface tension as 

P - Poo = 7V ■ n, (7) 
where the outward surface normal is given by 

n = . . ( - -^e, - -eJ , (8) 



and the divergence of the normal is 



V . n = — + -— (rn ) + -— ■ (9) 

dz r dr ^ r 86 



When linearized, this boundary condition becomes 



p'\r=R = -l{-^ + ^ + -^,^] ■ (10) 
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At the solid surface, the normal component of the velocity vanishes, and therefore 

u'-n = 0. (11) 

2.2. Normal modes 

Taking the divergence of the linearized Euler equation, Eq. [3l gives the Laplace equation 
for the pressure, 

VV = 0. (12) 
Considering normal modes for p', u' and ^' such that 

p'{r,e,z,t) =p{r,e)e''+''\ (13) 
u'{r,e,z,t) = u{r,e)e''+'''', (14) 
^'{e,z,t) =i{e)e''+^'\ (15) 

Eq. [12] becomes 

^ + 1^ + _ = (16) 

Substituting these normal modes into Eq. [3] relates the pressure to the velocity, 

dp I dp 

psvL = -TT^r T^ee - ikpe^, (17) 

or r o9 

allowing the boundary conditions to be expressed solely in terms of p. The boundary 
condition for the motion of the free surface, Eq. then becomes 

= (18) 

ps or 

The stress boundary condition on the free surface, Eq. [101 becomes 

Combining Eqs. [18] and [19] to eliminate ^ leads to 

s'p^4,(l-k"R'^i,)^- (20) 



pK^ \ 86'^ J dr 

Next, we need to express the no-penetration boundary condition at the wall in terms 
of p. Using the relations between p and u from Eq. [TTJ the no-penetration condition, 
Eq. [m may be written as 
dp I dp 

Tr'^r + -^ng = 0, (21) 
or r 06 

along the solid substrate. Finally, the condition that the contact points of the free 
surface are stationary requires 

1^0, (22) 

at the contact points. The complete eigenvalue problem to solve for the pressure field is 
therefore given by Eq. [161 together with the boundary conditions provided by Eqs. 
and I 
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We may solve Eq. [12] by the method of separation of variables, letting p{r, 9) = 
F{r-)G{6). Applying this separation gives 

F" F' _ G" _ ^ 

T^Vf'^-'^-^- ^^^^ 

We see that Eq. [23] for F is the modified Bessel differential equation, whose solutions 
are modified Bessel functions of the first and second kind, 

F{r) = C^h{kr) + C2K^{kr), (24) 

whereas the general solution to Eq. [23] for G is 

G{e) = Gs sin + G4 cos A^. (25) 

Applying the periodicity condition, that is, G{6) = G{6 + 2tx) to Eq. [25] gives that 
\n = n where n = 0,1, 2,.... Applying the condition that p is finite at r = to Eq. [2l] 
eliminates Kx as a solution for F. A general solution for p may then be written as 

00 

p (r, 6) = ^ In{kr) (A„ sin nO + B„ cos nO) , (26) 

n=0 

where An and -B„ are series of unknown constants to be determined by the boundary 
conditions. Applying the boundary condition at the free surface, Eq. [20], with r = R 
and 1^1 < 9c produces 

00 

^ IfiikR) {An sin nO + 5„, cos n9) (27) 

n=Q 

00 

= ^ V A; (1 - k^R'^ - I'(kR) (A^sinn^ + 5„ cosn^) . 

n=0 

The no-penetration condition at the lower wall, Eq. [2T], is given by 

where h is the vertical separation between the origin and the wall. Applying this 
condition using r = —h/ cos 9 and 9c < \9\ < n produces 



00 ^ 

n=0 ^ 



kh 



COS 9 



[An sin n9 — Bn cos n9) (29) 



n sin 6^ ^ / kh \ . 

— 7: [An COS n9 + Bn sm nt 

h \ cos 9 I 

n=0 ^ ^ 



At the contact points, r = R, \9\ = 9c, Eq. [22] becomes 

00 

^kl'nikR) [An sin n9 + Bn cos n9) = 0. (30) 



n=0 
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2.3. Eigenvalue problem 



Together, Eqs. [271 EH] and [30] represent a generalized eigenvalue problem, where the 
eigenvalue is the square of the growth rate, s^, and the eigenvector is composed of the 
series of constants An and Bn- Unlike many other eigenvalue problems which arise in 
separation of variables solutions, these equations cannot be solved in the usual manner, 
that is by multiplying by one of the modes, integrating over the domain of 6 and applying 
the orthogonality of the modes to generate analytical expressions for the constants An 
and Bn- This is because the boundary conditions in the physical problem are of mixed 
type. 

An approximate numerical solution to the eigenvalue problem may be sought by 
truncating the series at n = A^, and evaluating the boundary conditions at some number 
of discrete 9 = {9i, . . . , 6'2jv+i} to produce a linear system of 2N + 1 equations in 2N + 1 
unknowns {Aq is irrelevant). The set of 6 are picked to specifically include ±6c to ensure 
that the boundary conditions at those points are met. Formally, the eigenvalue problem 
of Eqs. [271— [30l may be rewritten as 



J2Anan{e) + Bnbn{9) = ^ A„a„ (e) + 5„/3„ (^) , 



n=0 



n=0 



where the functions a„, 6„, a„ and Pn are given by 



^rn{kR)k{l-eR^ -n^)smn9 \9\ < 9, 

krn(-^)smn9 + ^In(-^)cosn9 \9\ = 9, 

n V cos / n " V cos J II ' 

kl'nikR) smn9 



9r < \9\ < TT, 



^I'nikR) (1 - PR^ - n^) cosn9 



\9\ < 9. 



kl'i-^) cosn9 - ^U-^) smn9 \9\ = 9^, 

n\ cost)/ h "V cosOJ II 

kl'nikR) COS n9 9c<\9\<7r, 



(31) 



(32) 



(33) 



In{kR) sin n9 
a„ = < 




1^1 < 0,, 

9r < \9\ < TT, 



(34) 



and 



Pn 




\9\ < 9, 
\9\ = 9, 



< 



< IT. 



The approximate eigenvalue problem may then be written in matrix form as 



(35) 



ai(^i) 



«l(c'2Af+lJ 



0'n{Q2N+i) bo{92N+l) 



Ai 

An 
Bo 



B 



N 



(36) 
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Figure 4. Capillary instability of cylindrical segment on a hydrophilic stripe: Square 
of the dimensionless growth rate, / {'-^ / pR^) , as a function of the dimensionless 
wavenumber, kR, for various value of the apparent angle of the stripe, Oc = 90°, 
100°, 110°, 120°, 130° and 140°. We also reproduce as a dashed hue the growth rate 
for a free cylinder (classical Rayleigh-Plateau instability). 



ai ^1 



OCl{92N+l) 



OlN[tfl) 



2N+1} 



2N+1} 



(3N{d2N+l) 



A, 

An 
Bo 



B 



N 



The solutions to this generahzed eigenvalue problem may be found using Matlab or 
comparable software. 



3. Stability results 

The results of our linear stability calculation are illustrated in Fig. HJ where we plot the 
square of the dimensionless growth rate of the most unstable mode, s"^ / {•y / pR^) , as a 
function of the dimensionless wavenumber of the perturbation, kR, for different values 
of the apparent contact angle 6c (solid lines). We also plot for comparison the result for 
a free cylinder, i.e. 

2 7 mR) 



(1 - k'^R^) kR, 



(37) 



pR^IoikR) 

which is the classical Rayleigh-Plateau result (dashed line) An instability is 

possible only if there exists a value of k for which a mode of deformation satisfies > 0. 
We see in Fig. Hlthat, in accordance with previous work, the cylindrical segments with 

~ E a 



10 



11 



12 . The main result of the 



9c > 90° are always linearly unstable 
paper, as seen in Fig. HJ is the explicit calculation of the range of unstable wavenumbers 
(together with the associated growth rates) and in particular the result that this range 
reaches zero for 9c = 90°, which coincides with the limit of the stability domain. In 
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Figure 5. Normalized mode shapes, £,{0), as a fmiction of the angle 6 about the 
vertical direction, for a cylindrical segment with Oc ~ 110°, and wavenumber kR = 0.5. 
The three most unstable modes are plotted: s'^ / / pB?) = 0.008 (solid hne, blue; this 
is the only unstable mode), -2.539 (dashed line, green), -10.11 (dash-dotted line, red). 



other words, in the experiment of Ref . [5[ , as soon as the apparent contact angle reaches 
the critical value of 9c = 90° from below, the fluid becomes capillary unstable, but 
at a wavenumber that is close to zero, corresponding therefore to deformations with 
infinitely long wavelengths. Consequently, the experimentally observed bulges are the 
manifestation of a zero- wavenumber capillary instability [5]. We further plot in Fig. [5] 
the shape of the three most unstable modes for the position of the free surface ^{9) 
for 9c = 110° and kR = 0.5. Additional modes reflect higher spatial harmonics, all 
with negative values for s^. In no case does the value of become positive for any 
mode except the first, as is the case for the Rayleigh-Plateau problem [l|, Finally, 
we compare in Fig. [6] the results of our stability calculation (solid line) with that of the 
one-dimensional model of Ref. iJ] (dashed line) for 9c = 110°. We see that the one- 
dimensional approximation, although qualitatively similar to the the results of the full 
calculation, over-estimates both the growth rate and the most unstable wavenumber. 




Figure 6. Comparison between the inviscid stability calculation from the present 



paper (solid line) and the one-dimensional inviscid model from Schiaffino et al. 14 1 
(dashed line) for 9c — 110°. Square of the dimensionless growth rate of the most 
unstable mode, / {'-^ / pR^) , as a function of the dimensionless wavenumber, kR. 
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Figure 7. Fluid segment on a hydropliilic stripe on a wedge of opening angle 29c- (a) 
9c > 90° and (b) 9c < 90°. 



4. Cylindrical segment on a wedge 

Inspired by the results above, we consider now a different geometrical setup, as 
illustrated in Fig. [71 So far we have assumed the substrate to be horizontal. In that case, 
when the fluid volume is increased, both the apparent contact angle of the cylindrical 
segment on the surface and the cross-sectional shape of the fluid change. We now 
consider a setup where the contact angle is fixed, while the shape is allowed to change. 
Specifically, we consider a cylindrical segment of fluid on a hydrophilic stripe in a wedge- 
like geometry. The volume of the fluid and the opening angle of the wedge, 26c, are 
supposed to be such that the apparent contact angle of the fluid on the surface always 
remains to be 90°. As a result, the center of the cross section of the fluid segment is a 
circular wedge centered at the tip of the hydrophilic wedge (see Fig. [7]) . This geom etry 
is reminiscent of previous work on droplets in angular geometries, in both static 15|, [16 
and flowing [13] conditions. 

We consider therefore the capillary stability of the configuration illustrated in Fig. [7] 
with notation similar to the previous section. We now use a cylindrical coordinate system 
with the origin at the point of the wedge, and with 6 = aligned with one segment of 
the solid wall and 6 = 26c with the other. The boundary conditions in this geometry 
may be applied as follows. Applying the no-penetration boundary condition, Eq. [211 to 
the harmonic equation for the 6'-dependence of the solution, Eq. [251 Gives that C3 = 
and \n = wn I26c where n = 0, 1, 2, . . .. As before, requiring that p be finite at r = 
requires that C2 = 0. The form of the solution for p is then 

p(r, 6) = Y, Anh„ ( A;r ) cos ^ ■ (38) 



ra=0 



The boundary condition along the free surface, {r 
becomes 



R,0<6< 26c}, as given by Eq. 



S^J2^nhAkR) COS 



n=0 

00 

E 

n=0 



mrt 
~2dr 



(39) 



7 



kii- en' 



AJ'{kR) cos 



mT6 
~26^ 
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Figure 8. Capillary instability of cylindrical segment on a wedge: Square of the 
dimensionless growth rate of the most unstable mode, s^/ (7/pi?^), as a function of the 
dimensionless wavenumber, kR, for various value of the apparent angle of the stripe, 
0c 90°, 100°, 110°, 120°, 130° and 140°. We also reproduce in dashed line the growth 
rate for a free cylinder (classical Rayleigh-Plateau instability). 



At the contact points, {r = R, 9 = 0, 29c}, Eq. [22] becomes 

00 



(40) 



ra=0 



As above, we obtain a generalized eigenvalue problem for the eigenvalue and the 
eigenvector y4„. We can combine Eqs. HQ] and HQ] by writing them as 



J2 Anan{9) =S^Y1 '^nan{9), 
n=0 n=0 

where the coefficients a„ and a„ are defined as 



(41) 



ani9) 



Oin{9) 



7 



k{l- - i^) r^^{kR) COS ^ Q<9< 29c, 



h„{kR) cos ^ 0<9<29c, 



9 = 0oi 9 = 29r 



(42) 







or ^ = 29r. 



(43) 



As above, an approximate solution can be found by truncating the series aX n = N 
and evaluating Eq. Hi] at + 1 values of 9 which include ^ = 0. In matrix form, the 
truncated eigenvalue problem is now given by 



00(6*1) 



0-0(0 N+l) ■ ■ ■ ajv(^^jv+i) 

ao{9i) ■■■ aN{9i) 

OiQ{9N+l) ■ ■ ■ OiN{9N+l) 



A 



N 



(44) 



An 
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By numerically solving this eigenvalue system, we find that such a cylindrical wedge 
of fluid is stable as long as the opening of the wedge satisfies 6c < 90°. When 6^ > 90° we 
find unstable modes, and the square of their dimensionless growth rates are displayed in 
Fig. [8] as a function of the dimensionless wavenumber of the perturbation (solid lines). 
As above, we have included the classical Rayleigh-Plateau result (dashed line). We 
see that the growth rate of the most unstable modes in this setup (wedge-like stripe) 
are very similar to the ones obtained in the previous section (horizontal stripe). In 
particular, we recover the result that in the limit where the angle 9c — > 90°, the most 
unstable perturbation wavelength increases to infinity. 



5. Conclusion 



In this paper, we have studied the dynamics of the capillary instability discovered 
experimentally by Gau et al. In this work, it was shown that a circular segment 
of fluid located on a hydrophilic stripe on an otherwise hydrophobic substrate becomes 
unstable when its volume reaches that at which its apparent contact angle on the surface 
is ninety degrees. Instead of breaking up into droplets, the instability lead to the excess 
fluid collecting into a single bulge along each stripe. By performing a linear stability 
analysis of the capillary flow problem in the inviscid limit, we have first reproduced 
previous results showing that the cylindrical segment are linearly unstable if (and only 
if) their apparent contact angle is larger than ninety degrees. We have then calculated 
the growth rate of the instability as a function of the wavenumber of the perturbation, 
and shown that the most unstable wavenumber for the instability — the one which would 
therefore be observed in an experimental setting — decreases to zero when the apparent 
fluid contact angle reaches ninety degrees, allowing us to re-interpret the creation of 
bulges in the experiment as a zero- wawuumber capillary mstabilit; 0. A variafou of 
the stability calculation was also considered in the case of a hydrophilic stripe located on 
a wedge-like geometry. Since droplets of any shape can now be created experimentally 
using chemical substrate modification 18|] , the stability of more complex fluid topologies 
could be analyzed using a framework similar to the one developed here. 
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